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The electromagnetic scattering properties of topological insulator (TI) spheres are 
systematically studied in this paper. Unconventional backward scattering caused by 
the topological magneto-electric (TME) effect of TIs are found in both Rayleigh and 
Mie scattering regimes. This enhanced backward scattering can be achieved by intro¬ 
ducing an impedance-matched background which can suppress the bulk scattering. 
For the cross-polarized scattering coefficients, interesting antiresonances are found in 
the Mie scattering regime, wherein the cross-polarized electromagnetic fields induced 
by the TME effect are trapped inside TI spheres. In the Rayleigh limit, the quantized 
TME effect of TIs can be determined by measuring the electric-field components of 
scattered waves in the far field. 
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I. INTRODUCTION 


Topological insulators (TIs) are an emerging quantum phase in condensed matter 
physical - —. The key feature of TIs is that the electronic band structure resembles an 
ordinary insulator in the bulk, while gapless edge or surface states exist within the bulk 
energy gap! These edge states are protected by time-reversal symmetry. TI materials 
have been theoretically predicted and experimentally observed in various systems such as 
HgTe/CdTe and InAs/GaSb quantum well, and Bii_ x Sb x , Bi 2 Te 3 , Bi 2 Se 3 , Bi 2 Te 2 Sel'“. A 
novel quantized topological magnetoelectric effect (TME) is predicted in TIs that an applied 
electric field could induce parallel magnetization and an applied magnetic field could induce 
parallel electric polarization^. An additional term in the Lagrangian: AC = (@a/47r 2 )E • B, 
which gives rise to the “axion electrodynamics”, contains a complete description of the 
electromagnetic(EM) responses of Tls^d. Here a = e 2 /he is the fine structure constant, 
0 = (2 p + l)7r is the axion angle with p being an integer, E is electric field and B is mag¬ 
netic field. Owing to this additional Lagrangian, the constitutive relations for TIs should 
be modified as: D = eE — «B. H = B//z + cfE, where D and H are respectively the electric 
displacement and magnetic field strength, e and /j are respectively the permittivity and per¬ 
meability of TIs, a = Qa/n is proportional to the fine structure constant. These modified 
constitutive relations provide us a compact form to deal with the TME effect as the doping 
level is in the gap. Many unusual phenomena due to TME effect have been revealed^ - -’. 

The study of EM waves scattering by small particles and its application is an enormous 
field^. EM waves scattering by scatterers of TIs exhibits lots of interesting properties due to 
the TME effect. For instance, parity-violating light scattering under oblique incidence, and 
strong perturbation of dipole radiation near the TI spheres are demonstrated^. Broadband 
strong scattering in backward direction and interesting antiresonances are found for TI 
cylinders^. The antiresonance effect reveals that the TME effect is a surface effect essentially 
in the language of EM scattering. Moreover, the quantization of TME effect of TI cylinders 
can even be determined by measuring the electric-field components of scattered waves in the 
far fields at one or two scattering angles in the Rayleigh scattering limit—. 

In this article, we study light scattering by spheres of TIs. Based on Mie theory, we derive 
the scattering coefficients and scattering matrix for TI spheres analytically. Exotic strong 
backward scattering due to TME effect are also found in both Rayleigh and Mie scattering 
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regimes, similar to the case for TI cylinders. At certain frequencies, antiresonances of cross- 
polarized scattering coefficients of TI spheres are revealed, wherein the cross-polarized fields 
induced by TME effect are trapped inside TI spheres. In the Rayleigh limit, we propose 
a simple way to determine the quantized TME effect of TIs by measuring the electric-held 
components of scattered waves in the far held. 


II. MIE THEORY FOR TI SPHERES 



FIG. 1. Schematic view of the system under study. A TI sphere is placed at the origin. An EM 
plane wave is incident along z direction and polarized along x direction. 


The system under study is shown in Fig. [EJ We consider a TI sphere which is illuminated 
by a time-harmonic EM wave with frequency ui. The radius of the sphere is a. The dielectric 
permittivity and magnetic permeability for the TI sphere are (Ai,/ii), and those for the 
background medium are (£*,,//&)• This scattering problem can be solved analytically by the 
standard Mie theory. In the spherical coordinate system, (r, 9 , 0), the EM helds can be 
expanded by the vector spherical harmonics^: M^ n (fcr), M^ n (hr), N^(A:r) and N^ n (kr), 
where subscripts e and o denote respectively the even and odd modes with respect to the x 
axis, k is the corresponding wavevector ( k\ = y/einiui/c, kb = ^Eb^bu/c are the wavevectors 
in the TI sphere and background respectively), n is a integer, and / stands for which kind 
of spherical Bessel(Hankel) functions we use. We assume that the incident EM wave is 
Ej nc — E 0 e lkz e x , where E 0 is the amplitude. The scattered and internal electric helds can 
be expanded in terms of the vector spherical harmonics^: 
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( 1 ) 


E,c = 5] £„(-&„m 2>„ + S a„Ng'„ - 6”M® + «£'nS!.) 

n= 1 
n= oo 

Bw = X) E »( c » M iil. - «» N i‘iL + - «C nJ>) (2) 

n=l 


where E n = i n E 0 (2n+l)/n(n+l), {a n , &„} are the scattering coefficients that have one-to-one 
correspondence to the ones for the conventional dielectrics. However, two new scattering 
coefficients, {a^ 1 , bJJ}, or the so called cross-polarized, scattering coefficients arise from 
the TME effect. The superscripts 1 and 3 stand for the spherical Bessel function of the 


first kind and spherical Hankel function of the first kind, respectively. By matching the 


standard boundary conditions at r = a, and applying the modified constitutive relations 
D = eE — qB, H = B /fi'j- o-E inside the TI sphere, the scattering coefficients a n , b n , and 


the internal coefficients c n , d n are found to be: 

_ p b m 2 j n {rnx)[xj n (x)\(3 2 - pij n (x)[mxj n (mx)]' 
p b m 2 j n {mx ) [xh n (x)]'/3 2 ~ Pih n (x)[mxj n (mx)]' 
b _ Pijn(mx)[xj n (x)]' - p b j n (x)[mxjn(mx)]'/3 1 
Pijn(mx)[xh n (x)}' - p b h n (x)[mxj n (mx)}' 

_ Pljn(x)[xh n (x)}' ~ p 1 h n (x)[xj n (x)}' 

Pij n {mx)[xh n (x)]' - p b h n (x ) [mxj n (mx)]'/3 i 
d _ pimj n (x)[xhn(x)}' ~ pimh n (x)[xj n (x)\ 

p b m 2 j n (mx)[xh n (x)}'/3 2 - pih n (x)[mxj n (mx)]' 


(3) 

(4) 

(5) 

( 6 ) 


where j n and h n are respectively the spherical Bessel and Hankel functions of the first kind, 


x = k b a is the size parameter and m = yjsipi/y/e b p b is the relative refractive index; j3\ 
and f3 2 are auxilliary functions related to the axion angle, given by /3i — 1 + a 2 Xi and 
f3 2 = 1 — a 2 x 2 , where a = a^/pi/ei and 


Xi = 


X2 = 


p b m 2 j n (mx)[xh n (x)]' 


p b m 2 j n (mx)[xh n (x )]' - pih n {x)[mxj n {;mx)\ 
p b h n (x)[mxj n (mx)}' 


(7) 

( 8 ) 


pij n (mx)[xh n (x)}' - p b h n (x)[mxj n (mx)}' 

The cross-polarized scattering coefficients aTJ and b ^ 1 , corresponding to n-th order of electric 
and magnetic multipoles but with 90° polarization-rotation compared to conventional ones, 


are given by: 


\mxj n (mx)]' TI uTl _ 3n( mx ) T j 
m[xh n (x)]' n ’ n h n {x ) ° n 


(9) 
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where dJJ and c^ 1 are the cross-polarized internal coefficients of the TI sphere, given by: 

dn 1 = ic n , cl 1 = ax2d n ■ (10) 

Obviously, cross-polarized multipolar terms in the scattering fields: a^ 1 (electric) and b^ 1 
(magnetic) are caused directly by the internal multipolar terms: c n (magnetic) and d n 
(electric), respectively. This is exactly a strong manifestation of the TME effect that an 
applied electric (magnetic) field can induce magnetic (electric) field. For a topologically trivial 
insulator, the axion angle 0 = 0 , namely, <5 = 0 , therefore (3 1 = /3 2 = 1 , cl^ = bj/ = c^ 1 = 
d^ 1 = 0 , the scattering coefficients a n , b n and internal coefficients c n , d n are reduced to be 
the conventional ones^. 


III. SCATTERING COEFFICIENTS IN RAYLEIGH LIMIT 


We assume that Hi — Hb — 1 in this article since the materials considered here are non¬ 
magnetic. In the Rayleigh scattering limit (x <C 1 and nix -C 1), it can be shown that only 
the following scattering coefficients are in the order of x 3 : 


a\ = —i 


.2 3m 2 - 3 + 2a 2 /e b 3 5 

—zt^x 3 + 0(x 5 ) 


, .2 

Oi = i- 


3 3 m 2 + 6 + 2 a 2 /£b 
a 2 /e b 


3 3 m 2 + 6 + 2 a 1 /£b 


^ + 0(a; 5 ) 


af 1 = b\ 1 = i 


,ti 


‘lajy/Fb 


^ 27 — x 3 + 0(x 5 ) 


( 11 ) 

( 12 ) 

(13) 


3 m 2 + 6 + 2 a 2 /eb 

All other scattering coefficients are in the order of x 5 or higher, and can be hence neglected 
in the Rayleigh scattering limit. The scattering coefficients a\ and b\ correspond to the 
electric and magnetic dipolar terms, whereas af 7 and bj 1 correspond to the cross-polarized 
electric dipolar term and cross-polarized magnetic dipolar term respectively. 


IV. UNUSUAL BACKWARD SCATTERING 


Since the wavevector of the incident wave is along the z direction, the scattering direction 
e r and the forward direction e 2 define a scattering plane. The incident electric fields and 
the scattered electric fields in the far held (kr 1 ) can be related by a scattering matrix-: 


^11 

E s j_ 


i exp (ikr) 
kr 



Eq 

Ei± 


(14) 
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where E s y and E s j_ are respectively the components of scattered electric fields parallel and 
perpendicular to the scattering plane. The corresponding unit vectors parallel and perpen¬ 
dicular to the scattering plane are defined as: e s ||=e 0 and e s j_=-e,y The components of 
incident electric fields parallel and perpendicular to the scattering plane are respectively 
Ei\\ = Eq cos 0 and E i3 _ = E 0 sirup, which depend on azimuth angle. The elements of 
amplitude scattering matrix Sj (j= 1, 2, 3, 4) for a single TI sphere are given by the fol¬ 
lowing relations: Si = E“ i^TT) ( an7Tn + bnTn )> ^ = E“ ^TI) (' anTn + b n n n ), S 3 = 
En=r 5TI) ( a n Ir n ~ b l^n ), ^4 = ^T) {^n T n ~ «>»), where 7 T n = Pl/smd and 

T n = dP'/dO*. Pi are associated Legendre functions. For conventional dielectric spheres, 
the condition a ^ 1 = b ^ 1 = 0 leads to a diagonal scattering matrix, i.e., S 3 — S 4 , — 0. For 
TI spheres, however, scattering matrix is not a diagonal matrix in general since aJJ and bj/ 
may not vanish owing to the TME effect. 

The intensities of scattered waves, defined as energy flux per unit area, characterize 
scattering properties of the system. Based on the scattering matrix shown above, we’ll 
discuss the role that the dielectric background plays in the EM scattering. In the Rayleigh 
limit, only electric and magnetic dipolar terms dominate the EM scattering. The angle- 
dependent functions, 7 r n and r n for n = 1, are given by: 7 ly = l,Ti = cos 6 . Two different 
kinds of limit can be estimated in the following: 

(i) , ei — eb a. Therefore |ai| |af J | = \bj 1 \ |6i|, in other words, the contribution 

of the conventional electric dipole is much larger than those of the cross-polarized electric 
and magnetic dipoles induced by TME effect. The intensity of scattered waves is mainly 
determined by the diagonal terms of scattering matrix. Thus, the scattering patterns of TI 
spheres are similar to those of conventional dielectric spheres. 

(ii) , Ci — €b <C a. The refractive index of the background is very close to that of the TI 

sphere. In this case we have | aj 1 \ = |^P| |ai|, |6i|, namely, the conventional electric and 

magnetic dipoles are negligible compared to cross-polarized ones. We thus have: 


Is || - 
Is -L - 


( kr ) 2 

1 


|S 3 


9 1 


sm 0 = — 


a ^ 1 cos 6 — b{ 


TI I 2 


sm 


l^l 2 cos 2 0 = 


4 (kr ) 2 

9 1 1 . TI r\ TI |2 9 j 

b[ cos 6 — a( cos 0. 


(15) 

(16) 


(kr ) 2 4 (kr ) 2 

where we assume that the incident plane wave has an unit intensity, I s \\ and I s ± are intensities 
correspond to the polarized components with scattered electric-field vectors parallel and 
perpendicular to the scattering plane, respectively. The scattered intensity is contributed 
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FIG. 2. Distribution of scattered electric fields |E sca | 2 for TI spheres. The dielectric permittivity 
of the background medium is the same as that of TI spheres, (a) for Rayleigh scattering (x=0.01) 
and (b) for Mie scattering (x=l). The boundary of the sphere in (b) is indicated by a white circle. 
The incident wave propagates along z direction. 

mainly by off-diagonal terms of scattering matrix. Clearly, in the forward direction (9 = 
0), the scattered intensity reduce to zero for both I s m and I s ± since af 7 = 6 77 . In the 
backward direction (9 — n), however, the scattered intensity shows a strong backscattering 
stemming from the constructive interference of the cross-polarized electric and magnetic 
dipole radiation, i.e. the terms proportional to aj 1 and bj 1 . 

To show these unusual scattering patterns, we plot the distributions of scattered electric 
fields in Fig. [2j The dielectric permittivity of TI sphere and background medium are chosen 
to be £\ = £b = 30 . The resulting scattering patterns are completely different from those 
for the conventional media, showing a strong backward scattering for both Rayleigh and 
Mie scattering. It is noted that such enhanced or even complete backward scattering has 
also been found in other systems such as magnetio^ 1 ^, metallic^, or semiconducting^^ 
particles. However, unlike metallic spheres^ or semiconducting spheres with high refractive 
index—this enhanced backward scattering is a broadband effect for TI spheres since the 
TME effect is non-resonant intrinsically. 


V. ANTIRESONANCES 

In Fig. El the scattering and internal coefficients corresponding to dipolar modes as 
functions of size parameter x are shown. For scattering coefficients a\ and b \, there exist 
resonances and antiresonances. For instance, the resonances at x = 0.780 for |a,! | and 
x = 0.560 for |&i| correspond to the excitation of electric and magnetic dipolar modes 
respectively. Meanwhile, the internal coefficients di and C\ have similar resonant behaviors 
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FIG. 3. Scattering and internal coefficients as functions of the size parameter x with n = 1. 
The resonant peaks in (a) and (c) mainly come from the bulk scattering of TI spheres. The 
corresponding size parameters of antiresonances for cross-polarized coefficients are shown by dotted 
lines in (b) and (d). The dielectric permittivity of TI spheres is e = 30 and the axion angle is chosen 
to be 7 r. The background is vacuum. 

as a\ and b\ respectively, which indicates that the EM fields inside the sphere have also been 
strengthened at the corresponding frequencies. In contrast, at the antiresonances (x = 0.835 
for | ai | or x = 1.06 for |&i|) no such electric or magnetic dipolar radiation can be detected 
outside the sphere since the corresponding scattering coefficients vanish at these particular 
x. However, the electric and magnetic dipolar modes still exist inside the sphere, since the 
internal coefficients d\ and Ci are actually finite instead of zero as shown in Fig. [3fc). 

The scattering coefficients a n and b n and internal coefficients c n and d n for a TI sphere 
should be similar to those for a dielectric sphere with the same refractive index, since the 
hne structure constant is in the order of 10~ 2 . As a result, TME effect could not be mani¬ 
fested in the scattering coefficients a n , b n and internal coefficients c n: d n . However, it plays 
a key role in the cross-polarized coefficients a^ 1 , bJJ and c^ 1 , d^ 1 since these four coefficients 
should vanish in conventional dielectric spheres. Interestingly, there also exist resonances 
and antiresonances for a^ 1 and b^ 1 ■ The cross-polarized scattering coefficients vanish at the 
antiresonances. In order to guarantee a^ 7 = b^ 1 = 0, the condition j n (mx)[mxj n (mx)]' = 0 
should be satisfied. In other words, the antiresonances correspond to the roots of j n {mx ) = 0 
or [, mxj n (mx )]' = 0 exactly. In Fig. [3]^b), we take the dipolar mode as an example. The 



























FIG. 4. Field distributions in scattering plane <fi = 0 at antiresonances x = 0.501. In (a) and (b), 
the field components \Eg\ and | B^ of total fields contributed by the superposition of electric and 
magnetic dipolar modes are shown respectively. In (c) and (d), the components \E ( p\ and \Bg\ are 
trapped inside the sphere, which are contributed from the cross-polarized electric dipolar mode 
due to the TME effect, (e) and (f) show the field components \Eq\ and \B ( p\ of the electric dipolar 
mode arised from bulk scattering; while in (g) and (h), \Eg\ and |f?^| of the magnetic dipolar mode 
from bulk scattering are plotted. The field components Eg and B™ vanish at the boundary of the 
TI sphere. The physical parameters are the same as those in Fig. [31 The boundary of the TI 
sphere is shown by the dashed lines. 

first antiresonance at x — 0.501 corresponds to {mxji(mx)]' = 0, whereas the second an¬ 
tiresonance at x = 0.820 corresponds to ji(mx) = 0. Although the corresponding scattering 
fields are zero outside the sphere at the antiresonances, the fields inside the sphere are not 
necessarily zero. For instance, by checking the fields inside the TI sphere at x = 0.501, it 
is found that cross-polarized electric dipolar mode could be exited since dj 1 ^ 0 whereas 
cross-polarized magnetic dipolar mode does not exist since cf 1 = 0. On the contrary, for 
x = 0.820 we have cf 1 ^ 0 and dj 1 = 0 as shown in Fig. [3jd), which implies the existence of 
cross-polarized magnetic dipolar mode and the lack of cross-polarized electric dipolar mode. 

From the definition of the vector spherical harmonics^, the tangential components (i.e., 
eg and e^,) of N„ and M n are proportional to [mxj n (mx)}' and j n (mx) respectively. Thus, 
the underlying physics for the existence of anti-resonances of cross-polarized scattering co¬ 
efficients is that the tangential field components corresponding to N n or M n vanish at the 
boundary of sphere. Note that the electric(magnetic) fields of electric multipoles are con¬ 
tributed from N„(M n ), while those of magnetic multipoles are proportional to M„(N n ). 
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In other words, the corresponding tangential field components of the electric and magnetic 
multipoles should vanish at the boundary. In most cases, the electric and magnetic multi¬ 
polar fields yield from the bulk scattering are strong and the cross-polarized ones induced 
by TME effect cannot manifest themselves. However, TME effect can be clearly observed at 
two particular scattering planes, i.e., planes with 0 = 0° and 0 = 90°. In Fig. SJ we show the 
field distributions at scattering plane 0 = 0° for the first antiresonance at x — 0.501 which 
corresponds to a root of [mxji(mx)]' = 0. The upper panel of Fig. 0 ] shows the total fields 
distribution, the dominant field components are Eg and B<p , which are the superposition 
of the fields from electric dipole (polarized in x direction) and magnetic dipole (polarized 
in y direction) modes due to the bulk scattering. The field patterns are almost the same 
with those for conventional dielectric spheres. Interestingly, the field components £0 and 
Bg induced by cross-polarized electric dipole (polarized along y direction) are nonzero and 
trapped inside the TI sphere due to TME effect, whereas for conventional dielectric spheres 
these two components should be zero. As shown in Fig. |4](c) and (d), the cross-polarized 
fields induced by TME effect are completely trapped inside the sphere. Observation of the 
trapping effect can be made clearly since the incident and scattered waves don’t have these 
two components (£0, and Bg) at the scattering plane 0 = 0°. For conventional spheres^ 1 ^, 
the scattered waves are reduced to zero at antiresonances and thus the spheres become trans¬ 
parent, however, the corresponding field components of incident waves still exist inevitably. 

In contrast, at the antiresonances of TI cylinders, the electric multipolar modes (cor¬ 
responding to a n ) and magnetic multiple modes (corresponding to b n ) due to the bulk 
scattering can be excited separately via waves with different polarization, i.e., TE waves 
and TM waves respectively. As is well known, for TI spheres electric multipolar modes (a n ) 
and magnetic multipolar modes (b n ) will always be excited simultaneously. Therefore, at 
the antiresonances, the total fields of bulking scattering at the boundary can not be reduced 
to zero in general. However, the tangential components of the electric and magnetic fields 
for respectively the electric and magnetic multipolar modes, are zero as shown in Fig. He) 
and SKh) for the dipolar case with [rrixj\(rnx)]' = 0 satisfied. The vanishment of the tan¬ 
gential field components at antiresonances^ for electric and magnetic dipoles is an exotic 
phenomenon and also a strong evidence that the TME effect is intrinsically a surface effect. 
For the second antiresonance at x = 0.820 similar conclusions can be drawn. The field 
trapping for the cross-polarized magnetic dipole (note that c^ 1 ^ 0 here) inside the sphere 
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can also be observed. 


VI. DETERMINING THE FINE STRUCTURE CONSTANT IN 
RAYLEIGH LIMIT 


A quantitative determination of the TME effect can also be found in TI spheres similar 
to the case in TI cylinder—. We can conduct a Raylcigh-scattering experiment with incident 
plane waves as shown in Fig. [1] schematically. In the far field, both electric-field components 
of scattered waves |E S ||| and |E s j_| are measurable quantities and we can thus define a quan¬ 
tity R(9,(f>) = |E s i||/|E s _j_|. In the Rayleigh limit, the scattering is dominantly contributed 
from dipoles. Assuming that £& = Pt = Pi = 1 we have: 


R(M) = 


(a 2 — b cos 9) cos 0 — 3a(l — cos 6 ) sin 0 


(a 2 cos 9 — b) sin 0 — 3a (1 — cos 6 ) cos 0 


(17) 


where b = 3m 2 —3 + 2a 2 is basically a bulk parameter. To extract the fine structure constant, 


we firstly focus on the scattering angles that’re in the first quadrant, i.e. 9,(j) G (0,7 t/2 ). 
Obviously, a 2 is a small term compared to b, b cos 9 and hence could be neglected except in 


the vicinity 9 ~ tt/2. With these assumptions we have: 

_ bR(9, 0) sin0 — 6cos#cos0 . . 

3(1 — cos 9) [sin 0 — R(9, 0) cos 0] 

This offers a simple way to determine a if R(9, 0) is measured and the parameter m is known. 
Clearly, the parameter m is material-dependent and so does this measurement. However, if 
we set a special observation anlge (9, 0) = (90°, 0°), from Eq. (fl7l) we have: 


a = 3R(90°, 0°) 


(19) 


Owing to the complete suppression of bulk scattering in this particular angle, the to-be 
measured a now can be solely determined by the value of R(90°, 0°) which is material- 
independent. From the obtained a, it offers a way to measure the fine structure constant a 
since the axion angle is quantized. 


VII. CONCLUSIONS 

We’ve studied EM wave scattering by TI spheres systematically. The scattering coeffi¬ 
cients for TI spheres are derived analytically. We’ve found that the scattering patterns of TI 
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spheres showing strong backward scattering with impedance-matched background m ~ 1. 
The enhanced backscattering or the cancellation of the forward scattering is non-trivial since 
it is a strong manifestation of the TME effect where the EM scattering due to the surface 
of TIs play a key role. The strong backward scattering has a broadband feature and holds 
for both Mie and Rayleigh scattering regimes due to the non-resonant nature of TME effect. 
The anti-resonances in the Mie scattering regime are another kind of novel phenomenon of 
TI spheres, wherein the cross-polarized fields induced by TME effect are trapped inside the 
TI spheres. Finally, we proposed a simple way to determine the quantized TME effect of 
TIs in the Rayleigh limit by measuring the electric-field components of scattered waves in 
the far field. 
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